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A  PRISMATIC  FLOAT  HAVING  FREEDOM  IN  TRIM 
By  Robert  W.  Miller 


SUMMARY 


Equations  which  include  freedom  in  trim  are  derived  for  hydrodynamic 
impact  of  a  non-chine-immersed,  prismatic  float  forebody  having  a  V-bottom 
and  a  transverse  step.  These  equations  are  an  extension  of  previously 
published  fixed-trim  theory,  and  a  method  of  solution  is  indicated  by 
which  time  histories  of  vertical,  horizontal,  and  angular  displacement, 
velocity,  and  acceleration  can  be  obtained. 

Comparisons  of  specific  solutions  of  the  equations  with  corresponding 
fixed-trim  solutions  are  presented.  The  trends  and  deviations  noted  are 
similar  to  those  exhibited  by  a  like  comparison  of  experimental  data  for 
free-  and  fixed- trim  impacts. 


INTRODUCTION 


Previously  published  hydrodynamic  impact  theory  (references  1  to  3) 
has  been  based  on  the  concept  that  the  flow  about  an  immersing  seaplane 
float  or  hull  is  a  two-dimensional  phenomenon  occurring  in  transverse 
planes  fixed  in  space  and  oriented  normal  to  the  keel.  The  total  force 
on  the  float  is  obtained  by  summing  the  reactions  in  the  individual  flow 
planes  in  contact  with  the  float  and  applying  an  aspect-ratio  correction 
factor  to  account  for  the  end-flow  losses  which  exist  in  three-dimensional 
flow.  This  theory  has  made  use  of  the  simplifying  assumption  that  the 
trim  remains  fixed  throughout  the  impact.  Experimental  checks  of  this 
fixed-trim  theory  for  both  model  and  full-scale  hulls  have  been  presented 
in  numerous  reports  and  some  evaluation  of  the  empirical  factors  involved 
has  been  conducted. 

The  present  investigation  was  initiated  in  order  to  obtain  a  method 
for  determining  the  effect  of  freedom  in  trim  on  loads,  moments,  and 
motions  during  hydrodynamic  impacts  of  a  non-chine-immersed,  prismatic 
float  forebody  having  a  V-bottom  and  a  transverse,  step.  As  in  the  pre¬ 
vious  theoretical  presentations  the  wing  lift  was  assumed  to  be  equal 
to  the  model  weight  and  the  aerodynamic  moments  were  neglected  in  order 
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to  simplify  the  problem,  although  they  might  have  had  an  appreciable 
stabilizing  effect. 


The  present  paper  is  primarily  concerned  with  the  derivation  of 
equations  describing  the  vertical,  longitudinal,  and  angular  motions  and 
accelerations  of  the  'previously  described  float  during  hydrodynamic 
impact.  It  also  presents  comparisons  of  numerical  solutions  of  these 
equations  with  theoretical  fixed-trim  loads,  moments,  and  motions  under 
identical  initial  conditions,  and  a  comparison  is  made  of  experimental 
free-trim  and  fixed- trim  data. 


A 

a 


eg 


K 


l 

M 

eg 

Ms 


% 

•  • 

s 


s 

t 


V 

•  • 

X 


SYMBOLS 

hydrodynamic  aspect  ratio  (tan  p/tan  t) 
distance  between  center  of  gravity  and  step , 
distance  between • center  of  gravity  and  step, 
hydrodynamic  force 
moment  of  inertia  about  center  of  gravity  (pivot  point) 

constant  (|p^O)) 

wetted  keel  length 

moment  about  center  of  gravity  (bow- up  moment  is  positive) 

moment  about  step 

virtual  mass  of  float 

mss  of  fluid  in  a  flow  plane 

acceleration  of  center  of  gravity,  parallel  to  keel 
distance  of  flow  plane  forward  of  step,  parallel  to  keel 
time  after  contact 
velocity 

horizontal  acceleration  of  center  of  gravity 


parallel  to  keel 
normal  to  keel 
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Y 

y 

•  • 

z 

z 

;P 


horizontal  displacement  of  step  from  contact  point 
vertical  acceleration  of  center  of  gravity 
draft  at  step 

acceleration  of  center  of  gravity,  normal  to  keel 
penetration  of  float  into  a  flow  plane 
angle  of  dead  rise 


7  flight-path  angle  ^an"1  ^ 

P  mass  density  of  water 

t  trim,  radians ,  except  as  noted 

F(P)  dead-rise  function 

9(A)  aspect-ratio  correction  factor 

Subscripts: 

o  condition  at  time  of  water  contact 

n  direction  normal  to  keel 

p  direction  parallel  to  keel 


Where  units  are  not  specified  any  consistent  system  of  units  may  be 

used. 

Dots  over  x,  y,  z,  and  r  indicate  differentiation  with  respect 
to  time. 

METHOD  OF  ANALYSIS 


Physical  concepts.-  The  physical  concepts  on  which  this  analysis  is 
based  are  discussed  in  detail  in  references  1  and  3,  in  which  the  solu- 
ti°n  is  restricted  to  fixed  positive  trim.  Briefly,  the  theory  discussed 
m  these  references  was  based  on  the  concept  that  the  primary  flow  about 
an  immersing  seaplane  float  or  hull  occurs  in  transverse  planes  which 
were  considered  fixed  in  space  and  oriented  perpendicular  to  the  keel. 
Figure  1  is  a  sketch  of  a  prismatic  float  immersing  at  positive  trim. 
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A  given  flow  plane  is  indicated  in  figure  1(a)  at  a  distance  s  forward 
of  the  step  and  the  penetration  of  the  float  into  this  flow  plane  is 
represented  by  z.  The  motion  of  the  fluid  in  each  plane  was  treated 
as  a  two-dimensional  phenomenon  independent  of  the  other  flow  planes. 

The  total  force  on  the  float  is  obtained  by  integrating  the  reactions 
of  the  fluid  in  the  individual  flow  planes  in  contact  with  the  float. 

This  force  is  reduced  by  the  application  of  an  end-flow  (aspect-ra  10) 
correction  to  account  for  the  losses  which  exist  in  three-dimensional 
flow.  The  effects  of  buoyancy  and  viscosity  are  neglected,  since  in  an 
Impact  they  are  normally  small  in  comparison  with  the  inertia  forces. 

A  lift  force  eq.ua.l  to  the  weight  of  the  seaplane  is  also  assumed  to  act 

throughout  the  impact. 

These  concepts  are  followed  in  the  present  analysis  except  for  thfe 
modifications  necessary  to  introduce  freedom  in  trim.  In  order  to 
introduce  freedom  in  trim,  the  flow  planes  are  assumed  to  rotate  about 
instantaneous  centers  of  rotation  at  their  intersection  with  the  float 
keel.  Thus,  the  planes  are  permitted  to  maintain  their  orientation 
normal  to  the  keel  and  ah  the  same  time  to  have  no  translation  in  space 
along  the  line  of  the  keel.  . 

Equations  of  motion.-  In  this  analysis  both  the  float  being  con- 
sidered  and  the  forces  applied  to  it  are  symmetrical  about  a  vertical 
plane  through  the  keel  so  that  the  float  motions  can  be  resolved  as  a 
case  of  plane  motion.  The  equations  of  motion  for  the  float  (see 
fig.  1(a))  can  then  be  written  as 


Fp  =  -mfS  •  (-*-) 

Fn  =  -mpZ  (2) 


In  equations  (l)  and  (2)  the  quantities  S  and  Z  are  the  com-  , 
ponents  of  the  center-of-- gravity  acceleration  which  are,  respectively, 
parallel  and  perpendicular  to  the  keel.  These  accelerations  can  be 
expressed  in  terms  of  the  horizontal  and  vertical  components  of  the  step 
acceleration  and  the  center-of-gravity  accelerations  relative  to  the 
step  in  the  following  manner: 

s  =  X  cos  t  -  y  sin  t  -  ut2  .-  b*+  W 


A. 
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and 

Z  =  y  cos  t  +  x  sin  t  +  bT^  -  aT 


5 

(5) 


Since  the  viscous  forces  of  the  fluid  are  considered  small  enough 
to  he  neglected  in  this  analysis,  there  is  no  force  on  the  float  in  the 
direction  parallel  to  the  keel  and  hence  Fp  is  equal  to  zero.  There-  . 

fore,  from  equations  (l)  and  (4)  the  following  expression  can  he  obtained: 


-mf  (x  cos  t  -  y  sin  t 


(6) 


This  equation  can  then  be  solved  to  obtain  an  egression  for  x: 


•*  •  p  •• 

..  y  sin  r  ar  bT 

x  =  - +  -  +  - 

COS  T  COS  T  COS  T 


(7) 


Now  consider  equation  (2),  the  normal-force  equation.  The  normal 
force  on  the  float  is  the  sum  of  the  forces  exerted  by  the  individual 
flow  planes  reduced  by  an  end-flow- loss  correction  factor.  The  typical 
*  flow  plane  illustrated  in  figure  1(a)  is  at  a  distance  s  forward  of 

the  step  and  has  a  width  ds  parallel  to  the  keel.  The  float  has  made 
a  penetration  z  into  the  plane.  From  the  analysis  of  reference  3> 
the  virtual  mass  of  fluid  in  this  flow  plane  is 

mw  =  Kz2ds  (8) 

where 

K  =  |pnF(p)  (9) 

In  accordance  with  the  assumptions,  the  reaction  of  the  fluid  in 
the  plane  is  the  time  derivative  of  the  momentum  so  that 

dF  =  — (m^z)  =  — (Kz2ds  z)  =  K(z2£+  2zz2)ds  (10) 

*1  dt  dt 


* 
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The  normal  force  acting  on  the  float  can  now  he  obtained  by  inte¬ 
grating  over  the  wetted  length  and  multiplying  by  the  aspect-ratio  cor¬ 
rection  factor 


i=2  riy/sinT 

Fn  =  <P(A)  j  dF  =  Kqp(A)  /  (z^z  +  2zz^)ds 


(ID 


where  Fn  can  he  seen  to  he  a  function  of  z,  i,  and  z.  These  quan¬ 
tities  may  he  established  as  follows. 

It  can  he  seen  from  figure  1(a)  that  the  penetration  z  of  the 
float  into  a  given  flow  plane  can  he  written  as  a  function  of  the  vari¬ 
ables  y,  s,  and  r9  each  of  which  is  a  function  of  time.  Thus 


s  sin  t 


z  - 


COS  T  COS  T 


Differentiation  of  z  with  respect  to  time  gives 


(12) 


.  y  s  sin  t  yT  sin  t  st 

z  = - - +  - -p - — ^-5— 

COS  T  COS  T  COS^T  COS^T 


(13) 


and  a  further  differentiation  gives 


z  = 


s  sin  t  yV  sin  t  st  yi^  2yf^sin^T  2sT^sin  t 

-  -f  -  -  -  + - (-  - - -  -  - r -  + 


COS  T  COS  T  COS^T  COS^T  COS  T  COS^T 


COS 


3t 


2yt^sin  t  2s t 


cos^t 


p 

COS^T 
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The  derivatives  s  and  *s*  are  now  expressed  in  terms  of  x,  y, 
and  t  .  In  accordance  with  the  assumptions,  the  velocity  s  is  con¬ 
stant  over  the  wetted  length  and  is  simply  equal  to  the  negative  of  the 
float  tangential  velocity  Vpj  thus/  from  figure  1(b), 

1 

s  =  -Vp  =  y  sin  t  -  x  cos  r  (15) 

The  derivative  of  this  expression  is 


ds  ••  ••  •  •  •  • 

—  =  y  sin  t  -  x  cos  t  +  yr  cos  r  +  xt  sm  t 

dt 


Introducing  equation  (7)  gives 


s  =  yT  cos  t  +  xt  sin  t  -  ar^  -  hT 


(16) 


Equations  (12),  (13)/  and  (l4),  with  s  and  V  replaced  by  equa¬ 
tions  (15)  and  (l6),  respectively,  are  now  introduced  into  equation  (ll). 
Integration  is  then  performed  to  obtain  the  normal  force  in  terms  of 
motions  of  the  step.  Substituting  the  resulting  expression  for  normal 
force  into  equation  (2),  where  Z  in  this  equation  is  replaced  by  equa¬ 
tion  (5),  gives  the  following  equation  for  motion  of  the  center  of  gravity 
normal  to  the  keel: 


Kcp(A) 


12  sin^T  COS^T 


•y1#*  sin3T  cos  t  -  sin  r  cos  t)  +  t^(16  sin^T  - 


6  sin^T  +  2)J  +  y3 1^(4  sin^T  cos2t)  +  V(4b  sin3T  cos^t)  + 


t^(4a  sin3T  cos2t)  +  xt(20  sin^T  cos^t)  4-  yT(l2  sin.T  cos3t 


20  sin  t 


c°s5t)J  +  y^ Jl2  sin^T  cos^T(y  cos  t  +  x  sin  t)^J 

f  y  bT  sin  t  .0  ..  at^sin  t\ 

mf  (— - +  - +  bT  -  aT  +  - -  J 

Vpos  T  COS  T  COS  T  / 


(17) 
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In  order  to  solve  the  third  equation  of  motion,  the  total  moment 
exerted  by  the  individual  flow  planes  about  the  center  of  gravity  must 
be  determined.  To  make  an  approximate  correction  for  end-flow  losses 
this  moment  is  multiplied  by  the  aspect-ratio  factor  9(A),  where  the 
assumption  is  made  that  this  over-all  correction  factor  can  be  applied 
to  the  moment  as  well  as  to  the  normal  force.  Equation  (3)  thus  be 
written  as 


Mcg  =  9(A)  / .  (s'  -  a)dF  =  IcgT 


which,  with  the  use  of  equation  (ll),  becomes 


MCg  =  Kqp(A) 


y/  sin 


(z2z  +  2zz2)s  ds  -  aFj 


Proceeding  as  with  the  normal-force  equation  gives  the  following  equation 
for  angular  motion  about  the  center  of  gravity: 


K9(A) 

60  sin^T  cos^ 


-«y 

t  L  — 


Q  *  O 

(5  sin-^T  cos  t  -  2  sin  t  cos  t)  +  t^(25  sin 


19  sin2T  +  6)J  +  y*|v*(5  sin2T  cos2t)  +  r(5a  sin2T  cos  t  - 

20a  sin^T  cos  t  +  5b  sin3r  cos2t)  +  sin^T  -  85a  sin^T  - 

10a  sin  t)  +  xt(25  sin2T  cos2t  -  35  sin2T  cos^t)  +  £4(15  sin  t  ,cos3 

35  sin  t  cos^tTJ  +  y3|jy(20a  sin3r  cos2t)  -  *-r(20ab  sinV  cos2t)  - 


T2(20a2sini|'T  cos2t)  _  xT(l00a  sin^T  cos2t)  +  yr(A0a  sin2T  cos3t  - 

100a  sin^T  cos3T)  +  20  sin2T  cos^T(y  cos  t  +  x  sin  t)^J  - 
y2^60a  sin3r  cos^T^y  cos  r  +  x  sin  T)qj|’  =  IcgT 
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Equations  (7),  (17),  and  (20)  are  the  equations  of  motion  for  the 
float  and  involve  the  three  unJsnowns  x,  y,  and  t.  These  equations 
can  he  solved  numerically  for  the  time  histories  of  the  three  variables  > 
and  their  derivatives;  see  appendix  A  for  a  step-by-step  iteration  pro¬ 
cess.  "  Involved  in  the  solution  of  these  equations  are  the  dead-rise 
function  (see  equation  (9))  and  the  aspect-ratio  correction  factor  (see 
equation  (ll)),  which  are  discussed  in  the  subsequent  section. 

Factors  F((3)  and  cp(A).-  In  order  to  perform  the  numerical  solu¬ 
tions  expressions  must  be  obtained  for  the  dead-rise  function  F(p)  and 
the  aspect-ratio  correction  factor  cp(A)  contained  -in  the  normal-force 
and  moment  equations. 

An  expression  for  the  dead-rise  function  F(|3),  which  gives  the 
variation  of  the  effective  two-dimensional  fluid  mass  with  dead  rise, 
was  given  by  Wagner  in  reference  4  as 

m  =  - 1)2 

which  is  equivalent  to  (f(pj]^  of  reference  3. 

Although  this  relationship  has  not  heen  experimentally  verified 
for  impacts  of  floats  having  low  dead-rise  angles,  it  has  been  found 
to  be  in  substantial  agreement  with  test  data  obtained  with  floats 

lO  Q 

of  22g-  to  40  angle  of  dead  rise  (see  reference  3)  • 


On  the  basis  of  experiments  with  vibrating  plates  in  water  (refer¬ 
ence  5) ,  Pabst  derived  the  following  aspect-ratio  expression  for  approxi¬ 
mating  the  three-dimensional  virtual  mass  from  the  virtual  mass  computed 
on  a  two-dimensional  basis: 


9(A) 


where  A  is  the  aspect  ratio  of  the  equivalent  vibrating  plate.  If 
this  reduction  in  the  virtual  mass  is  assumed  to  be  determined  by  the 
shape  of  the  intersected  area  in  the  plane  of  the  water  surface,  then 
the  application  of  Pabst* s  data  to  V-bot tom  floats  results  in  the 
expression 


9(A) 


tan  t 

1 - - - 

2  tan  P 
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THEORETICAL  RESULTS 


Reduction  to  fixed,  trim.-  The  free-trim  equations  were  derived 
from  the  same  basic  assumptions  as  the  fixed-trim  equations;  they  must, 
therefore,  be  reducible  to  the  fixed-trim  equations.  This  reduction  is 
accomplished  in  appendix  B  by  considering  the  trim  to  be  a  constant 
instead  of  a  variable,  so  that  the  terms  containing  t  and  t  drop 
out.  It  is  then  demonstrated  that  equations  (17)  and  (20)  reduce  to  a 
form  identical  with  the  equations  for  normal  force  and  step  moment  in 
references  1  and  6,  respectively.  Equation  (7)  is  also  shown  to  reduce 
to  a  simple  fixed-trim  relation. 

Comparison  of  free-  and  fixed- trim  theoretical  solutions.-  In  order 
to  indicate  the  effect  of  freedom  in  trim  both  free-  and  fixed-trim 
theoretical  solutions  were  made  for  a  typical  hull  for  three  different 
sets  of  landing  conditions.  The  dimensions  and  inertia  values  of  the 
full-scale  float  considered  are  given  in  table  I.  The  theoretical 
initial  conditions  for  the  three  impact  cases  treated  are  listed  in 
table  II.  (The  experimental  values  also  listed  in  table  II  are  discussed 
in  a  subsequent  section. ) 

Figure  2  presents  time-history  comparisons  of  solutions  of  the 
free-trim  equations  with  fixed-trim  solutions  for  the  three  cases  of 
theoretical  initial  conditions  given  in  table  II.  The- quantities 
plotted  are  the  vertical  and  angular  displacements,  velocities,  and 
accelerations. 

The  free-trim  solutions  were  obtained  by  the  method  described  in 
appendix  A  and  the  fixed-trim  solutions  were  obtained  by  the  method  of 
reference  7.  The  so-called  ¥  curves  for  fixed  trim  shown  in  the  plots 
were  obtained  by  using  equation  (3)  but  with  the  fixed-trim  total  moment 
from  equation  (B7)  instead  of  the  free-trim  moment.  These  r  curves 
can  be  interpreted  as  a  measure  of  the  applied  moment  since,  from  equa- 

tion  (3),  V  =  -C-£  and  Icg  is  the  constant  value  for  the  free-trim 
Icg 

float. 


EXPERIMENTAL  RESULTS 


To  demonstrate  the  effects  of  freedom  in  trim  experimentally,  both 
free-  and  fixed-trim  tests  of  a  float  model  were  made  in  the  Langley 
impact  basin.  The  measured  motions  are  compared  in  this  section. 
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No  direct  comparison  “between  the  experimental  and  theoretical 
results  can  be  made  because  of  differences  between  the  float  used  in 
the  tests  and  that  assumed  in  the' theory.  The  float  used  in  the  tests 
had  a  pointed  step ,  flared  chines,  and  an  afterbody;  whereas  the  theory 
applies  to  a  float  forebody  having  a  transverse  step  and  an  infinitely 
wide  V-bottom. ‘  Since  these  float  tests  represent  the  only  free-trim 
impact-basin  data  now  available,  however,  the  results  are  included 
herein  for  the  qualitative  comparisons  that  can  be  made. 


Apparatus  and  instrumentation.-  The  Langley  impact  basin  and  testing 
procedure  are  described  in  reference  8.  Pertinent  dimensions  of  the 
float  model  used  in  the  experiments  are  given  in  table  I.  This  model, 

except  for  the  differences  outlined  previously,  is  a  -  -scale  model  of  • 

2 

the  float  considered  in  the  theoretical  solutions. 


In  the  free-trim  tests  the  float  model  was  attached  to  the  dropping 
linkage  of  the  launching  carriage  at  two  main  pivot  points;  these  sup¬ 
ports  were  on  a  transverse  line  which  passed  through  the  point  corre¬ 
sponding  to  the  center  of  gravity  of  the  airplane  from  which  the  float 
was  patterned.  For  the  fixed-trim  tests  a  third  support  point  was 
located  about  20  inches  aft  of  the  main  supports  on  the  longitudinal 
center  line  of  the  float. 

For  the  free -trim  tests  the  model  was  supported  only  at  the  two 
main  points  and  was  held  at  fixed  trim  until  just  prior  to  contact  by 
means  of  a  locking  mechanism.  After  contact  the  model  was  free  to 
rotate  in  pitch  about  the  main  supports  over  a  trim  range  of  -6.5° 
to  22.5°.  Beyond  those  limits  angular  displacement  of  the  float  was 
restrained  by  two  shock  struts  which  were  coupled  to  the  float  by  means 
of  telescoping  tubes,  one  60  inches  forward  and  one  60  inches  hft  of 
the  main  pivots.  The  buffer  action  of  the  shock  struts  extended  the 
trim  range  about  5°  in  each  direction  before  a  stop  was  reached. 

Two  strain-gage  accelerometers  of  the  same  type  of  construction 
were  electrically  connected  to  obtain  angular  acceleration  directly. 
These  accelerometers  were  located  on  a  longitudinal  line  passing  through 
the  axis  of  rotation  and  at  a  distance  of  6  feet  forward  and  6  feet  aft 
of  this  axis.  A  control-position  transmitter  was  adapted  to  the  equip¬ 
ment  to  measure  angular  displacement. 

A  standard  NACA  three -component  accelerometer  was  used  to  obtain 
the  vertical  component  of  acceleration  of  the  float.  It  had  a  natural 
frequency  of  21  cycles  per  second  and  a  critical  damping  of  0.8.  Other¬ 
wise,  the  standard  instrumentation  as  described  in  reference  8  was  used. 

Comparison  of  free-  and  fixed- trim  experimental  results.-  Figure  3 
presents  the  comparison  of  free-trim  and  fixed-trim  experimental  results 
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for  cases  I  and  II  listed  in  table  II.  Again  vertical  and  angular 
displacements,  velocities,  and  accelerations  are  compared;  the  values 
shown  apply  "to  full-scale  conditions • 

It  cap  be  noted  in  table  II  that  the  initial  vertical  velocity  of 
the  fixed-trim  run  of  case  I  was  higher  than  that  of  the  corresponding 
free- trim  run.  The  time  histories  for  the  fixed- trim  run  were  there- 
fore  adjusted  as  follows  to  give  more  comparable  results  in  figure  i^a;,. 
Since  the  initial  flight-path  angles  of  the  two  runs  were  not  grea  y  ^ 
different  and  the  initial  trims  were  the  same,  the  velocity  time  history 
of  the  fixed- trim  run  was  scaled  by  the  ratio  of  the  free-trim  ni  1a 
velocity  to  the  fixed-trim  initial  velocity,  the  acceleration  time-, 
history  was  scaled  by  the  square  of  this  ratio,  and  the  time  was  scaled 
by  the  reciprocal  of  the  ratio.  No  corresponding  adjustments  were  made 
for  case  II. 


DISCUSSION 


Initial  conditions.-  Table  II  presents  the  initial  conditions  for 
■both  the  free-  and  fixed-trim  theoretical  solutions  and  also  for  the 

free-  and  fixed-trim  test  runs.  The  test  runs  were  made  with  a  --  scale 

model  but  the  initial  conditions  are  given  in  terms  of  full-scale  values 


The  initial  conditions  for  the  theoretical  solutions  of  cases  I 
and  II  were  chosen  to  correspond  to  the  free-trim  experimental  runs  o. 
those  cases  to  make  the  results  as  comparable  as  possible.  The  experi¬ 
mental  fixed-trim  runs  were  chosen,  from  among  those  available,  to 
correspond  as  closely  as  possible!  to  the  respective  free-trim  runs. 

The  initial  conditions  for  case  III  represent  a  more  severe  impact 
such  as  might  occur  in  a  second  or  third  contact  during  a  landing  run 
or  in  a  contact  against  the  flank  of  an  advancing  wave. 


Case  I  -  The  first  set  of  conditions  treated  (case  i)  represents 
an  impact  at  moderate  initial  values  of  flight-path  angle  and  trim.  In 
both  the  theoretical  (fig.  2(a))  and  the  experimental  (fig.  3(a))  results 
the  vertical  motions,  trim,  and  angular  velocity  for  the  free-trim  con¬ 
dition  differ  only  slightly  from  the  corresponding  motions  in  the  fixed- 
trim  condition.  The  only  angular  acceleration  or  moment  of  any  size  is 
exhibited  by  the  experimental  free-trim  run  (t  ~  0.15),  "but  this  momen 
is  associated  with  the  immersion  of  the  flared  chines.  Thus,  for  the 
initial  conditions  of  case  I  no  significant  difference  appears  to  exist 
'between  fixed-  and  free-trim  impacts. 
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Case  II,-  Case  II  has  approximately  the  same  initial  flight-path 
angle  as  case  I  hut  has  a  higher  initial  trim.  The  vertical  motions 
for  the  free-trim  condition  are  again  observed  to  approximate  those 
for  the  fixed-trim  condition  in  both  the  theoretical  (fig.  2(b))  and 
experimental  (fig.  3(t>))  curves.  This  impact,  however,  resulted  in  an 
appreciable  change  of  attitude  and  a  moderate  nose-down  moment. 

The  initial  trends  of  the  theoretical  and  experimental  angular 
motions  are  similar.  However,  the  Immersion  of  the  flared  chines  prob¬ 
ably  caused  large  loads  on  the  float  forward  of  the  center  of  rotation 
and  therefore  a  large,  nose-up  moment  as  indicated  by  the  positive  peak 
of  the  experimental  t  curve.  These  loads  would  result  in  large 
deviations  of  the  subsequent  motions  of  the  test  float  from  those  of 
the  float  assumed  in  the  theory. 

Case  III.-  Case  III  represents  a  high-flight-path-angle,  low- trim 
condition  such  as  might  occur  in  an  impact  on  the  flank  of  a  wave.  The 
curves  represent  the  "effective"  values  referred  to  the  inclined  plane 
of  the  water  surface  as  in  reference  9  and  are  analyzed  as  an  extreme 
condition  of  smooth -water  impact. 

In  this  case,  the  deviations  between  the  fixed-trim  and  free-trim 
results  are  much  greater  than  those  previously  described;  freedom  in 
trim  in  fact  reduces  appreciably  the  maximum  values  of  y  and  t.  The 
large  values  of  negative  vertical  velocity  of  the  free -trim  solution 
compared  with  the  fixed-trim  solution  late  in  the  impact  appear  to  be 
consistent  with  the  large  values  of  trim  attained. 

Probably  the  most  serious  effect  observed  for  the  free-trim  solu¬ 
tion  is  the  high  trim  and  large  positive  value  of  angular  velocity  at 
the  end  of  the  impact.  This  could  result  in  a  stalled  condition  between 
impacts  or  could  lead  to  extreme  initial  conditions  for  a  subsequent 
impact.  However,  these  conditions  should  be  somewhat  restricted  by  the 
presence  of  an  afterbody  and  by  aerodynamic  moments  not  taken  into 
account  by  the  present  theory. 

No  impact-basin  data  have  been  obtained  for  conditions  as  severe 
as  those  in  this  case.  An  indication  of  the  motions  to  be  expected, 
however,  can  be  obtained  from  an  impact  of  a  four-engine  flying  boat, 
the  data  for  which  have  not  been  published.  The  initial  conditions  for 
this  impact  were  substantially  the  same  as  those  of  case  III  and  the 
resulting  motions  were  also  very  similar.  Of  particular  interest  in 
this  seaplane  impact  is  the  increase  in  trim  from  a  small  initial  value 
(about  3°)  to  about  15°  with  positive  angular  velocity  at  the  time  of 
exit  from  the  water.  The  acceleration  records  contained  large  oscilla¬ 
tions,  due  probably  to  structural  vibrations,  which  rendered  them 
unsuitable  for  direct  comparison  with  the  theory;  however,  the  peak 
values  of  the  faired  curves  for  both  y  and  t  agree  roughly  with  the 
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corresponding  values  for  the  free-trim  solution  of  case  III.  Thus, 
although  direct  comparison  is  not  practicable,  this  flying-boat  impact 
indicates  that  the  initial  conditions  for  case  III  are  within  the 
practical  range  and  that  the  theoretical  results  which  are  obtained 
for  this  case  are  at  least  consistent  with  flight  experience. 

General  discussion.-  In  general,  for  cases  I  and  II,  the  theoreti¬ 
cal  free-trim  vertical  motions  approximate  the  fixed-trim  vertical 
motions  to  about  the  same  degree  as  was  found  experimentally.  Also, 
the  trends  of  the  theoretical  angular  motions  are  similar  to  those  of 
the  experimental  data  except  where  deviation  is  to  be  expected  because 
of  differences  in  body  shape.  Thus,  for  impacts  having  moderate  ^initial 
conditions  the  theoretical  and  experimental  free- trim  results  appear 
to  be  in  agreement. 

For  the  more  extreme  initial  conditions  of  case  III,  there  is,  of 
course,  more  deviation  of  all  variables  from  the  fixed-trim  solution. 

The  deviations  which  occur  are,  however,  consistent  among  themselves 
and  with  general  experience.  Also,  the  reductions  of  maximum  applied 
load  and  moment  which  occur  are  to  be  expected. 


CONCLUDING  REMARKS 


Equations  of  motion  are  derived  and  presented  for  free-trim  hydro- 
dynamic  impact  of  a  V-bottom,  transverse  step,  prismatic  float  forebody. 
The  equations  are  an  extension  of  previously  published  fixed-trim  hydro- 
dynamic  impact  theory  and  are  chiefly  based  on  the  same  concepts  and 
assumptions.  A  method  of  solution  of  the  equations  is  also  presented 
which  gives  the  results  as  time  histories  of  vertical,  horizontal,  and 
angular  displacement,  velocity,  and  acceleration. 

The  free-trim  equations  are  shown  to  be  reducible  to  the  fixed- 
trim  case.  Moreover,  comparisons  of  specific  solutions  of  the  free- 
trim  equations  with  fixed-trim  solutions  for  the  same  initial  conditions 
and  also  with  some  experimental  data  have  shown  that  the  computed 
results  are  reasonably.  The  free-trim  solutions  and  experimental  data 
exhibit  similar  trends  and  deviations  from  the  fixed-trim  case,  and 
changes  of  attitude  and  applied  moment  are  consistent  with  previous 
experience. 

Although  these  facts  do  not  completely  validate  the  equations, 
they  at  least  indicate  that  reasonable  and  consistent  results  can  be 
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obtained  by  their  use.  Further  experimental  investigation  should 
provide  a  better  evaluation  of  the  method. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  February  13,  1952 
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APPENDIX  A 

METHOD  OF  SOLUTION 

Before  the  numerical  solution  can  he  performed,  equations  (17) 
and  (20)  must  he  transformed  into  forms  better  suited  to  the  method  of 
solution  used.  Equation  (7)  is  already  in  the  form  required  for 
solution. 

Equation  (17)  can  he  rewritten  as 

At  +  By  =  C  (Al) 

and  equation  (20)  as 

Dr  +  Ey  =  F  (A2) 

where 


Kcp(A)  [y^(4  sin^T  cos  t  -  sin  t  cos  t)  +  y^(kb  sin^T  cos^t)J 

12  sin^T  cos^t 
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C  =  - 


S' 

- . . . <  y^-|T2(l6  sinV  -  6  sin2T  +  2)J  + 

12  sin3r  cos5r  _ 


■3jT^(4a  sin^T  cos2t)  +  xt(20  sinV  cos^t)  4*  yr(l2  sin  t  cos^t  - 


20  sin  t  cos^t)  +  y2^  sin^r  cos^T(y  cos  t  +  x  sin  t 


(aT  sin  t  . p 

- - +  bT 

.  COS  T  J 


Kcp(A) 

60  sin^T  cos5 


£5(5  sin^T  cos  t  -  2  sin  t  cos  t)  + 


y^(5a  sin2T  cos  t  -  20a  sin^T  cos  r  +  5b  sin^T  cos2t)  _ 
y3(20at  sin^r  cos^t^J  -  Icg 


K<p(A) 

60  sin^T  cos^t 


Jy^(5  sin2T  cos2t)  -  y3(20a  sin^T  cos2t)J  (A7) 
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F  = 


Kcp(A) 


60  sin^T  cos^r 


•5^(25  sin^T  -  19  sin^T  +  6)J 


^[^(35a  sin^T  -  85a  sin^T  -  10a  sin  t)  +  xt(25  sin2r  cos2t 
Ln2T  cos^t)  +  yr(l5  sin  t  cos^t  -  35  sin  t  cos^ioj  - 


35  sin 


•3 ( 20a2  s in^  t  cos2t)  +  xT(l00a  sin^T  cos2t)  -  yT(40a  sin2r  cos^t 


100a  sin 


cos^t)  -  20  sin2 r  cos^^y  cos  t  +  x  sin  t)^]  - 


y^j&Ja  sin^T  cos^y  cos  t  +  x  sin  t) 


(A8) 


Equations  (Al)  and  (A2)  can  Toe  solved  simultaneously  to  give 


..  AF  -  DC 

y  =  ' - 

■  AE  -  DB 


and 


(A9) 


••  EC  -  BF 

T  =  - 

EA  -  BD 


(A10) 


which,  together  with  equation  (7), 

•  •  •  2  •  • 

..  y  sin  t  aT  Tor 

x  =  -  +  - - +  - — 

COS  t  COS  T  COS  T 


are  the  equations  to  he  used  in  the  numerical  solution. 
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The  numerical  solution  can  he  performed  by  extrapolation  of  the 
plots  of  y  and  r  against  time  by  any  desired  function  to  obtain  an 
approximate  value ^  at  the  end  of  the  next  time  interval.  The  other  quan¬ 
tities  (x,  x,  x,  y,  y ,  t ,  and  t)  are  then  computed  and  the  process 
iterated.  By  this  method,  time-history  plots  of  each  of  the  nine  vari¬ 
ables  are  obtained.  The  present  solutions  were  made,  however,  by  the  , 
o 

Kutta  ^  method  (reference  10 )  which  is  more  readily  adaptable  to  the 

8 

automatic  computing  machine  used. 
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Similarly,  the  right-hand  side  of  equation  (20)  is  equal  to  the 
moment  about  the  center  of  gravity  as  in' equation  (3)  and,  if  the  trim 
derivatives  are  set  equal  to  zero  in  the  left-hand  side,  the  equation 
becomes 


Kcp(A) 


60  sin^T 


cos- 


-  5: 

)T  i— 


5y  t 


2 

sin  t 


2 

COS  T 


-  20ay3y  sin^r  cos^t  + 


(SOay^sin^T  costly  cos  t  +  x  sin 


(b6) 


which  reduces  to 


•M 


Kq)(A)y^y  Kcp(A)y3(y  cos  t  +  x  sin  t)2  Kcp(A)ay3*y 


Cg  2  -3  ' 

12  sin  t  .cos-'t 

Kcp(A)ay2(y  cos  t  +  x  sin  t)^ 


3  sin^T  cos  t 


3  sin  t  cosjt 


(B7) 


Sin  T  COS  T 


After  equations  (B3)  and  (B4)  are  substituted  into  equation  (B7),  the 
moment  becomes 


Mcg  = 


Kcp(A) 

ySr 

3  sin^T  cos^t 

T" 

+  y3cos^T(y  cos  t  +  x  sin  t)^ 


K»(A)y3  ^ 

O 

3  sin  t  cos  t 


Kcp(A)y%n2 
sin  t  cos  t 


(b8) 


and,  hy  introducing  equation  (B5)  and  equations  (9)  and  (11) -of  refer¬ 
ence  3,  equation  (B8)  becomes 
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Mcg  +  aJ?n  = 


Kcp(A) 

k*  • 

y  T 

3  sin^r  cos^t 

_k 

+  y3(y  +  Ki  cos  t) 


(B9) 


where  is  .  K  in  reference  3. 

From  equations  (3^)  'and.  (26)  of  reference  6  the  moment  ahout  the 
step  is  i 

Ms  =  Mcg  +  aFn  (BIO) 


Substituting  equation  (BIO)  into  equation  (B9)  gives  the  moment  ahout 
the  step 


Ms 


Kcp(A) 

k»  • 

y  t 

3  sin^T  cos^t 

_4 

y3(y  +  K]_  cos  t) 


(Bll) 


which  is  identical  with  equation  (16)  of  reference  6. 

\ 

Equation  (7)  can  also  he  reduced  in  this  manner  to  give 

x  =  y  tan  t  (B12) 

In  the  fixed-trim  condition  the  step  motions  and  center-of-gravity 
motions  are  the  same.  Equation  (B12)  is,  therefore,  equivalent  to 

/ 

X  •=  Y  tan  t  (B13) 


Equation  (B13)  can  readily  he  deduced  from  figure  l(c). 
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PHYSICAL  CHARACTERISTICS  OF  FULL-SCALE  HULL  AND  g  -  s6aLE 

»  •  1 

EXPERIMENTAL  MODEL 


Full- scale 
Lull 


Experimental 

model 


scale 


Beam,  ft  . .  . . .  *  *  *  * 

Forebody  length,  from  centroid  of  step,  ft  .  . 

Over-all  length,  ft  . .  • 

a  (from  centroid  of  step),  ft  . 

b  (normal  to  keel),  ft  •  ••••••••• 

Weight,  lb . . . 

mj,  lb-sec2/ft  . .  . 

Icg,  lb-sec2-ft 

P,  deg  .  .  .  . . 


6.33 
19.11 
35-67 
2.28 
6. 06 
14,000 

434.78 

19,4io 

20 


3.17 

9.56 

17.83 

l.l4 

3.03 

1,750 

54.35 

607 

20 


LExperimental  values  changed  to  full-scale 
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(a)  Sketch  of  float  during  impact. 


(h)  Velocity  components.  (c)  Components  of  center-of-gravity 

acceleration. 


Figure  1.-  Schematic  representation  of  impact. 
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Figure  3*-  Concluded. 
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